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Comment on "Variation of the superconducting transition temperature of hole-doped 
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We point out the incorrect derivation of the gap equation in X.-J. Chen and H. Q. Lin [Phys. 
Rev. B 69, 104518 (2994)] within the interlayer tunneling (ILT) model for multilayered cuprates. 
There, the local structure in k-space of the ILT effective interaction has not been taken into due 
account when the ILT model is generalized to the case of n layers per unit cell. This is a specific 
characteristic of the ILT model that, apart from giving rise to a highly nontrivial k-dependence 
of the gap function, is known to enhance the critical temperature T c in a natural way. As a 
consequence, we argue that Chen and Lin's results cannot be employed, in their present form, for 
a quantitative interpretation of the high-pressure dependence of T c in Bi-2212, as is done by X.-J. 
Chen et al. cond-mat/0408587 to appear in Phys. Rev. B]. Moreover, when the generalization of 
X.-J. Chen et al. cond-mat/0408587 is applied to the case n = 2, it fails to reproduce the original 
ILT gap equation. However, a more careful analysis of the ILT model for multilayered cuprates, 
taking into account the nonuniform hole distribution among inequivalent layers, has been earlier 
suggested to describe the observed pressure dependence of T c in homologous series of high-T c 
cuprates. 
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PACS numbers: 74.62.-c, 74.72.-h, 74.62. Fj, 74.20.-z 



In Ref. [I], Chen and Lin reconsider the dependence of 
T c on doping and on the number of layers in a homol- 
ogous series of multilayered high-Xj, cuprates within the 
interlayer tunneling (ILT) model However, in deriv- 
ing their gap equation, Chen and Lin erroneusly neglect 
the intrinsic local structure in momentum (k) space of 
the effective ILT coupling. This is a specific character- 
istic of the ILT model, which is known to give rise to 
highly nontrivial features in the k-dependence of the gap 
function already for a bilayer complex Moreover, 
a local term in the gap equation has been shown to pro- 
vide a lower bound for T c at all dopings, which is the 
precise way in which the ILT mechanism enhances T c Q . 
The consequences of such an incorrect analysis of the 
ILT model arc both qualitative and quantitative. There- 
fore, the recent use of Chen and Lin's results to interpret 
the high-pressure dependence of T c in Bi-2212 p| can 
be questioned. In this context, we point out that a more 
careful analysis of the ILT model for layered cuprates has 
been presented elsewhere @, and successfully applied to 
study the pressure dependence of T c in homologous series 
of layered cuprates, by explicitly taking into account the 
inhomogeneous hole-doping in inequivalent layers 0, ■ 

Superconductivity in the high-T c layered cuprates is 
characterized by (i) a non-monotonic dependence of T c 
on the overall hole-doping 5; (ii) a monotonic increase 



of T c with the number of layers n, for moderately low 
n (n < 3). While (i) is a generic consequence of the 
quasi-bidimcnsional nature of these compounds (see e.g. 
Ref.0), 

the latter fact has suggested that coherent tun- 
neling of superconducting pairs between adjacent Cu02 
layers may considerably enhance T c 0. Within the 
ILT model, it is postulated that strong in-plane corre- 
lations forbid coherent hopping of single particles be- 
tween adjacent CuC>2 planes. Such a restriction is re- 
moved when accessing the superconducting state, where 
interlayer Josephson tunneling of Cooper pairs is allowed. 
This results in a net gain in kinetic energy, as compared 
to the normal state. Thus, within the ILT model, su- 
perconductivity is stabilized via a kinetic mechanism, as 
opposed to conventional BCS superconductivity, where 
the enhancement in kinetic energy is overcompensated 
by a reduction in the potential energy [Tjj . 

However, after its original formulation more than a 
decade ago Q, the relevance of the ILT mechanism at 
least for single-layer cu prat es has been called into ques- 
tion by experiments [llLll2j . Recently, Chakravarty et al. 
[l3j have revived the ILT model in connection with multi- 
layered cuprates. There, ILT needs not be the sole source 
of superconducting condensation energy. Charge carri- 
ers require a 'seed' in-plane interaction to form Cooper 
pairs in a given symmetry channel, before they can ac- 
tually tunnel between adjacent layers Such in-plane 
interaction would then provide the missing condensation 
energy [l3j . 
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Moreover, it has been suggested that the competition 
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with a 'hidden' order parameter, such as a d-dcnsity-wave 
(dDW) , could be responsible for the downturn of T c 
with n, for n > 3. Indeed, in multilayered cuprates, due 
to the different proximity to the 'charge reservoir' blocks, 
experiments 1 1 5| as well as density functional theory cal- 
culations [lfij revealed a nonuniform hole-content distri- 
bution between inner and outer layers. Since this usually 
places inner (outer) layers in the undcrdoped (overdoped) 
region of the cupratc phase diagram, competition with 
the dDW order would be stronger in inner layers than in 
outer layers, thus depressing T c with increasing n. Hy- 
drostatic pressure could then be used to tune both the 
overall hole-content content and its distribution among 
inequivalent layers, thus inducing an 'exchange of roles' 
between inner and outer layers with respect to the onset 
of superconductivity 0, which is observed as 'kinks' in 
the pressure dependence of T c in layered cuprates . 

The effective Hamiltonian considered by Chen and Lin 
(see also 00) is 
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where £k is the in-plane quasiparticle dispersion mea- 
sured with respect to the chemical potential // (// = /i 
for all layers, in Ref. 0), and c£ a is a quasiparticle cre- 
ation operator with wave-vector k and spin a on layer t. 
It should be emphasized that in Eq. the first interac- 
tion term (Vkk') pertains to a single layer and governs the 
overall symmetry of the order parameter (i.e., ci- wave, if 
Vkk' = Vg k g k >, with g k = ^{cosk x -cosk y )) H3, while 
the second term applies to adjacent layers (££'), and is 
local in k-space, with Tj(k) = j^Tj(cosk x — cosfc y ) 4 0. 
This enforces momentum conservation for the intcrlaycr 
pair tunneling process. (The effect of k-spacc broaden- 
ing of the ILT kernel, e.g. due to impurities, has been 
considered in Ref. ITsh 

A straightforward mean-field analysis of Eq. Q for a 
bilayer complex and an in-plane superconducting insta- 
bility in the d-w&ve channel yields the gap equation 0: 



A k = 



Aog k 



i-Tj(k) Xk ' 

where Ao is determined self-consistently from 
V 
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Here, Xk = (2-Ek) 1 tanh(/3£'k/2) is the pair suscepti- 
bility at the inverse temperature (3 = (fefiJ 1 ) \ = 
VSk + |Ak| 2 is the upper branch of the superconducting 
spectrum, and N is the number of lattice sites. 

Eqs. {2J an d © should be immediately compared and 
contrasted with Eq. (9) in Ref. (for a multilayered 
complex) and Eq. (1) in Ref. (for a bilayer complex). 



Even without going into the subtleties of the more gen- 
eral derivation for an n-layered complex (for which, see 
Refs. 0, or with the competition among several in- 
plane pairing channels (Ref. |3|), it is apparent that the 
gap function within the ILT model, Eq. is character- 
ized by a local prefactor [1 — Tj (k)xk] _1 which, albeit 
linked self-consistently to Ao via Eq. J3J), is responsible 
of most of the quantitative and qualitative features of 
the model. Such a structure is missing in Refs. 
Even though the actual symmetry of the gap function 
is independent of the ILT kernel, and is rather deter- 
mined by the d-wave nature of the in-plane coupling, the 
ILT mechanism endows the gap function with a nontriv- 
ial structure in k-space 0, which has been shown to be 
consistent with ARPES results 0. Moreover, the 'renor- 
mahzed' pair susceptibility in the summand of Eq. (ffi . 
viz. Xk l— * Xk/[1 — Tj(k)xk], which is due to the local 
ILT tunneling amplitude, gives rise to additional, alge- 
braic divergences in the energy dependence of the inte- 
grated pair susceptibility, as opposed to the logarithmic 
one, typical of BCS theory This is directly respon- 
sible of the enhancement of T c within the ILT model. In 
particular, in the case of a bilayer complex, one analyti- 
cally finds a lower bound for T c as 
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where nearest (t) and next-nearest (t') hopping have been 
assumed, and fi± = —At + 4t', fij- = At + At' , and /^vh = 
—At' denote the bottom, the top of the band, and the 
location of the Van Hove sing ularity, respectively 0. 

On the contrary, the ILT kernel enters Chen et al.'s 
gap equation in Eq. (9) of Ref. and Eq. (1) of Ref. 
as an additional contribution to the non-local in-planc 
coupling term, i.e. it amounts to defining another in- 
plane interaction, with no reference to interlayer tunnel- 
ing. This same observation applies to the general case of 
an n-layered complex. In that case, the gap equation for 
each layer should also contain a local contribution due 
to the ILT mechanism between adjacent layers (again, 
absent in Ref. 0, with an ILT renormahzed pair suscep- 
tibility x k /[! - T j( k )xL] for each layer 00, with 
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which can be further simplified in the limit of uniform 
hole-content in all layers (as is tacitly assumed in Ref.[l|). 
In analogy to the bilayer case, the condition 



min[l - Tj(k)x k ] = 

k 



(0) 



then implicitly defines a lower bound T* 1 for the critical 
temperature corresponding to the onset of superconduc- 
tivity in the given layer I. Therefore, for nonuniform 
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hole-content among inequivalent layers, as is the case for 
the multilayered cuprates |l5l Il6l |. one can estimate a 
lower bound to T c as max^T* £ . A nonuniform distribu- 
tion of the overall hole-content among inequivalent layers 
can be conveniently described by means of appropriate 
models Q . This then enables us to identify whether the 
superconducting instability first sets in in inner or outer 
layers. One finds a crossover as function of the overall 
hole-content Q, which has been related to the observed 
'kinks' in the pressure dependence of T c in several layered 
cuprates Q- 

In conclusion, we have pointed out an incorrect deriva- 
tion of the gap equation(s) for layered cuprates within the 
ILT model pi 0| for the general case of n superconduct- 
ing layers per unit cell. This in turn leads to a failure 
to capture most of the qualitative and quantitative fea- 



tures of the theory, both for bilayered and multilayered 
compounds. As a consequence, the theoretical analysis 
of the high pressure data in Ref. is not consistent with 
the ILT mechanism. On the other hand, a more careful 
analysis of the ILT model when taking into account 
a nonuniform hole-content distribution among inequiva- 
lent layers Q, is indeed able to reproduce the observed 
pressure dependence of T c in multilayered cuprates Q . 
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